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AN IMPROVED METHOD OF COLLOCATION FOR THE STRESS ANALYSIS 


OF CRACKED PLATES WITH VARIOUS SHAPED BOUNDARIES* 

By J. C. Newman, Jr. 

Langley Research Center 

SUMMARY 

An Improved method of boundary collocation was developed and applied to the two- 
dimensional stress analysis of cracks emanating from, or in the vicinity of, holes or 
boundaries of various shapes. The solutions, presented in terms of the stress -intensity 
factor, were based on the complex variable method of Muskhelishvili and a modified 
boundary-collocation method. The complex-series stress functions developed for simply 
and multiply connected regions containing cracks were constructed so that the boundary 
conditions on the crack surfaces are satisfied exactly. The conditions on the other 
boundaries were satisfied approximately by the modified collocation method. This 
improved method gave more rapid numerical convergence than other collocation tech- 
niques investigated. 


INTRODUCTION 

In the life of a structure subjected to cyclic loads, cracks may initiate at and propa- 
gate from geometric discontinuities (holes, cut-outs, edges, or flaws). In designing 
structures for fatigue and fracture resistance the stresses around these sites must be 
calculated if unexpected failures are to be avoided. In recent years the use of high- 
strength (crack sensitive) materials for greater structural efficiency has resulted in a 
series of aircraft service failures due to the presence of small cracks. The rates at 
which these cracks propagate and the size of the crack that causes failure are strongly 
influenced by the shape of the discontinuity and the type of loads applied near it. 

The crack -tip stress-intensity factor (restricted to small-scale yielding) can 
account for the influence of component configuration and loading on fatigue-crack growth 


The information presented herein is based in part upon a thesis entitled "Stress 
Analysis of Simply and Multiply Connected Regions Containing Cracks by the Method of 
Boundary Collocation," offered in partial fulfillment of the requirements for the degree of 
kfe.ster of Science, Virginia Polytechnic Institute, Blacksburg, Virginia, May 1969. 
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and static strength (refs. 1 to 4). Many investigators (refs. 5 to 12) have obtained theo- 
retical stress-intensity factors for cracks growing near stress concentrations or 
boundaries. However, for cracks in complicated structures, stress-intensity factors can- 
not always be obtained by closed-form analytical methods. They are often obtained by a 
series solution (such as boundary collocation) which, to be useful, should converge rapidly 
to the proper answer. 

The objective of this paper is to present an improved method of collocation for 
calculating the influence of various boundary shapes on the stress -intensity factor in two- 
dimensional linear elastic bodies and to apply this method to a number of boundary -value 
problems involving cracks. 

Stress-intensity factors were calculated for cracks emanating from a circular hole 
in an infinite plate and cracks emanating from a circular hole in a finite plate, a crack 
near two circular holes in an infinite plate, a crack emanating from an elliptical hole in 
an infinite plate, and a crack in a finite plate. These configurations were subjected to 
either internally or externally applied loads. 

SYMBOLS 
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coefficients in the complex stress functions 
half crack length 

half axis of the elliptical hole perpendicular to the plane of the crack 

coefficients for jth pole in the complex stress functions 

half axis of the elliptical hole, in the plane of the crack 

distance from the centerline of the crack to the center of the circular hole 

stress-intensity correction factor 

resultant force per unit thickness acting in the x and y directions, 
respectively 

sth influence function for nth coefficients 
resultant forces and displacements 
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H 

i = rr 

K 

Kt 

M 


N 

P 

P 

R 

r 

S 

u,v 

W 

z 


a 

/3 


specified resultant forces or displacements on the boundary 
height of rectangular plate 

stress-intensity factor 
stress-concentration factor 

number of points at which the error in the boundary conditions was 
evaluated 

number of coefficients in each series stress function 
concentrated force per unit thickness acting in the y direction 
pressure 

radius of the circular hole 

minimum radius of curvature for the ellipse 

stress on the external boundary 

displacement in the x and y direction, respectively 

width of rectangular plate 

complex variable, z = x + iy 

location of pole in x-y plane, Zj = Xj + iy^ 

ar^le between the x-axis and the normal to a boundary 

constant in equation (1) 

coordinate along the contour of a boundary 



K 


material constant 


A. ratio of applied stresses on separate portions of a boundary 

/i Lame’s constant (shear modulus) 

u Poisson’s ratio 

(Tj^ normal stress at a boundary 

Ox normal stress in x direction 

CTy normal stress in y direction 

shear stress at a boundary 
Txy shear stress 

complex stress functions 

Subscripts: 
j,m,n,q,s indices 

METHOD OF BOUNDARY COLLOCATION 

Boundary collocation is a numerical method used to evaluate the unknown coeffi- 
cients in a series stress function, such as those developed in appendix A. The method 
begins with a general series solution to the governir^ linear partial differential equation. 
Certain terms are eliminated from the series by conditions of symmetry. The series is 
then truncated to a specified number of terms. The coefficients are determined through 
satisfaction of prescribed boundary conditions. The series solution finally obtained 
satisfies the governing equation in the interior of the region exactly and one or more of 
the boundary conditions approximately. 

Various techniques have been used to satisfy boundary conditions. Conway (ref. 13) 
determined the coefficients from the criterion that the boundary conditions be satisfied 
exactly at a specified number of points on the boundary. Hulbert (ref. 14) and Hooke 
(ref. 15) selected the coefficients so that the sum of the squares of the stress residuals 
was a minimum for a specified number of points on the boundary. These two techniques 
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have been used to analyze the stress state around a crack in a rectangular plate (refs. 11 
and 14). Bowie (ref. 12) used a modified mapping-collocation technique to calculate the 
stress-intensity factors for a crack in a circular disk. 

The present approach combines the complex variable method of Muskhelishvili 
(ref. 16) with a modified boundary -collocation method. The modification requires that 
the resultant forces on the boundary be specified (in a least-square sense) in contrast to 
previous work in which the boundary stresses were specified. (See appendix B.) 

The resultant forces and displacements are expressed in terms of the complex 
stress functions ^(z) and J2(z) (see appendix A) as 


/3j^(z)dz + Jn(z)dz + (z - z)4>{z) = f(x,y) + ig(x,y) 


( 1 ) 


The bars denote the complex conjt^ates. For resultant forces (/3 = 1) acting over the 
arc ? - ?Q on the boundary 



For convenience, the location of Cq was selected as the intersection of the boundary 
with either the x- or y-axis. For displacements (/3 = -k) at a point ? on the bovindary 


2|n(u + iv) = -|l(x,y) + ig(x,y)^ ^ 


(3) 


For the case of plane strain k = 3 - 4v and for plane stress /c = 


3 - V 
1 + 


In general, the truncated complex stress functions can be expressed as 


N 

<p(z) = ^ An0j^(z) 
n=0 

> (4) 

N 

S2(z) = ^ BnOn(z) 

n=0 J 

where and Ojj are power-series functions of z. From equations (1) and (4) the 
expressions for f and g are 
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( 5 ) 


where Fgj^ and Gg^ (s = 1,2) are influence functions derived from and 
for the nth coefficients. 


Having developed the expressions for f and g in terms of the unknown coeffi- 
cients An and Bn, it is necessary to select a method for evaluating the coefficients 
from the specified boundary conditions. The two methods developed for this purpose 
were (1) direct specification of boundary forces and/or displacements and (2) least- 
squares specification of boundary forces and/or displacements. 


Direct Specification of Boundary Forces and/or Displacements 

In specifying the resultant forces and/or displacements on the boundary, equa- 
tions (5) were evaluated at N equally spaced points on the boundary and the resulting 
equations were solved on a computer for the unknown coefficients. Further details on 
the computer and matrix solutions are given in the section entitled "Application of the 
Boundary -Collocation Method." 


Least-Squares Specification of Boundary Forces and/or Displacements 

In general, the computed boundary values are in error at any given point 
because the series stress functions were truncated. The square of this error is 
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where f^ and gQ are specified boundary values at point The coefficients are 

then evaluated by minimizing the squares of the errors at a specified number of 
points M on the boundary: 
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= 0 
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= 0 


( 7 ) 


J 


where q = 1,2,. . .,N. Equations (7) result in a set of 2N linear algebraic equations 
for the unknown coefficients An and Bji 




N N 

^ ^nq'^n "*■ ^ ^q®n ~ ^q 
n=0 n=0 




N N 

^ ^qn^n + ^ '^nq®n = 
n=0 n=0 


(8) 


where 


M 

®nq “ ^ (^In^lq ^In^lq) 
m=l 


^q 


M 

2 (Wlq + G2nGlq) 

1 ' ill 
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M 

Z (^2nF2q + G2nG2q) 

m=l “ 

M 

X (^oFlq + goGlq)^ 
m=l 

and 

M 

"^q= Z (^oF2q + goG2q) 

m=l 

The second partial derivative of the square of the error with respect to the unknown 
coefficients was positive, indicating a definite minimum. 

APPLICATION OF THE BOUNDARY-COLLOCATION METHOD 

The improved collocation method was used to analyze a number of crack problems. 
The configurations investigated were grouped into four categories: (1) cracks emanating 

from a circular hole in an infinite plate, (2) a crack near two circular holes in an infinite 
plate, (3) cracks emanating from an elliptical hole in an infinite plate, and (4) a crack in 
a finite plate. In each category several boimdary conditions were investigated. 

In applying the method, it was necessary to specify the points on the boundary at 
which the error equation was to be evaluated. The locations on the circular and elliptic 
boundaries were specified by dividing the angle 0 q into M equal increments (A0) 

(see fig. 1). This procedure was selected because more points are concentx’ated on sec- 
tions of the boundary which have smaller radii of curvature and, hence, larger stress 
gradients. For simplicity, this procedure was also employed on the external boundaries. 
In general, the value of M used in the solution of the various boundary -value problems 
was twice the total number of unknown coefficients in the series stress functions. For 
the case of cracks emanating from the elliptical hole in an infinite plate and from the 
circular hole in a finite plate, 160 coefficients were used in the stress functions. In 
all other cases 90 coefficients were used (except where noted). 

A digital computer was used to solve the resulting equations using either single 
precision (14 significant digits) or double precision (29 significant digits). The equations 
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were solved by a routine which employed Jordan’s method (ref. 17). Using single preci- 
sion, the computer time was 1 and 4 minutes for solving 90 and 180 equations, respec- 
tively. Double precision required approximately twice as much time as sii^le precision. 
Double precision was used only in the cases involving the finite plate. 

The results are presented in terms of a correction factor F which accounts for 
the influence of the various boundaries on the stress-intensity factor for a crack in an 
infinite plate. The correction factor is defined as the ratio of the stress-intensity factor 
for the particular case to that for a crack in an infinite plate subjected to the same 
loading. A table of numerical values of F and the stress -intensity equation are given 
in Tables I to VII. 

Cracks Emanating from a Circular Hole in an Infinite Plate 

Remote biax ial stress .- For the case of cracks emanating from a circular hole, 
four collocation techniques were used to obtain the unknown coefficients in the stress 
functions. (See eqs. (All) of appendix A.) The conditions at infinity were satisfied by 
adding to equations (All) the n = 0 terms in equations (AlO). The results of each of 
these techniques are compared in figure 2. The curves were obtained by specifying 
either stresses at discrete points on the boundary (stress equations) or resultant forces 
along the bovindary (force equations), or by using a least-squares procedure with the 
stress equations or force equations. For the least-squares procedure employing resul- 
tant forces, two convergence curves are presented. One is for the case in which the 
number of points M was twice the total number of coefficients in the stress functions 
(2N is the total number of coefficients in the stress functions); the other is for the case 
in which the number of points considered was five times the total number of coefficients. 
Convergence is seen to be rapid in both cases. For the least-squares procedure 
employing stress equations, the number of points considered was twice the total number 
of coefficients. All techniques seemed to converge to the same value as the number of 
coefficients increased. The value of F was approximately 2 percent higher than 
Bowie's approximate solution (ref. 5) for the same configuration and loading. The 
boundary stresses also converged to their specified values as the number of coefficients 
increased. However, the analysis using the least-squares procedure with the force 
equations converged more rapidly than the other techniques. 

The correction factors for three states of biaxial stress are shown in figure 3. 

The solid curves show the results obtained in the present investigation. The circles are 
Bowie's solution and were obtained from a table given in reference 18. The agreement is 
considered good. 
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A solution for other values of X was constructed by superposition of the two 
cases, X = 0 and X = 1. The stress-intensity equation is 

K = - X)Fq + XF^ (9) 

where Fq and Fj are the correction factors for X = 0 and X = 1, respectively. 

The calculations for X = -1 confirmed equation (9) to greater than eight significant 
digits. However, equation (9) is only valid for positive values of stress intensity. Nega- 
tive stress intensities would indicate crack closure and the problem would become one of 
calculating contact stresses. 

Internal pressure .- In a configuration like that in the preceding section, internal 
pressure p was applied to the hole boundary and Xp to the crack surfaces. The 
stress functions for this case are given by equations (All). 

The correction factors for X = 0 and X = 1 are shown in figure 4. For the case 
in which no pressure is applied to the crack surfaces (X = 0) the correction factor 
modifies the stress-intensity factor for wedge-force loading on the crack surfaces 
(ref. 18). Here the wedge-force loading 2pR is that produced by the internal pressure. 
The other case modifies the solution for a uniformly pressurized crack. A solution for 
other values was constructed by superposition of these two cases. The stress -intensity 
equation is 


j^^2£R(i_x)f +Xp\^Fj (10) 

V jra 

where Fq and Fj^ are the correction factors for X = 0 and X = 1, respectively. 
Equation (10) is also valid only for positive values of stress intensity. The results indi- 
cate that the influence of the hole on the stress-intensity equation may be neglected for 
values of a/R greater than 2.5. 

A Crack Near Two Circular Holes in an Infinite Plate 

Crack located between the holes subjected to remote uniaxial stress.- In the case 
of two circular holes the stress functions, equations (A13), were used with poles located 
at Zj = ±id (see appendix A). The conditions at infinity were satisfied by adding to 
equations (A13) the n = 0 terms in equations (AlO). 

The correction factors are shown in figure 5 for several values of d/R. At small 
d/R values the holes had a pronounced effect on the correction factors. However, for 
d/R values greater than 10 or a/R values greater than 3 the influence of the holes may 
be neglected. 
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In the limit as the crack length approaches zero, the correction factor is equivalent 
to the local-stress-concentration factor at the origin. In terms of the remote 

stress and the local stress the stress-intensity factors are 

K = S/mF = Oy/lfa (11) 

from which the correction factor is given as the ratio of Oy to S. The circle plotted 
on the ordinate axis is the local-stress-concentration factor at the origin as given in 
reference 19 for a value of 4- = 2. The agreement is considered good. 

ix 

Cracks approaching the holes subjected to remote uniaxial stress .- In this case the 
stress functions, equations (A 13), were used with poles located at z. = ±d. Again, the 
conditions at infinity were satisfied by adding to equations (A13) the n = 0 terms in 
equations (AlO). 

The correction factors are shown in figure 6 for several values of d/R. The 

factors are plotted against the nondimensional crack length ^ ■ . The correction 

a - R 

factors increase from their initial value at a = 0 to large values as the crack 
approaches the edge of the hole. Again, the correction factor at a = 0 is equal to the 
local-stress-concentration factor at the point midway between the two holes. 

As the crack intersects the boundary of the holes, the concept of the stress- 
intensity factor K no longer exists. However, the stress-concentration factor K-p at 
the edge of the hole (x = R + d) does describe the severity of the notch. The stress- 
concentration factors are shown in figure 7 as a function of R/d ratios. The stress - 
concentration factors were compared with those calculated from an "equivalent" 
elliptical-hole solution as shown by the dashed line. The elliptical hole had a minimum 
radius of curvature equal to the radius of the circular holes and had an overall length 
equal to 2(R + d). The stress-concentration factors for the elliptical hole were only 
slightly lower than those calculated for the case of two circular holes connected by a 
crack. Thus, the elliptical -hole solution may be used to approximate this configuration. 

Cracks Emanating from an Elliptical Hole in an Infinite Plate 

Remote uniaxial stress .- For the case of cracks emanating from an elliptical hole 

(see fig. 8) the stress functions, equations (A13), contained multiple poles Zj on either 

the X- or y-axis. The n = 0 terms in equations (AlO) were added to equations (A13) in 

|_ 

order to satisfy the boimdary conditions at infinity. For the case of — = 0.5 and 2, the 

b ^ 

value of J in equations (A13) was 4 and for — = 0.25 and 4 the value was 16. These 
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values were determined by trial and error. The poles were always located on the major 
axis of the ellipse at the origin, at the center of the minimum curvature, and equally 
spaced between these points. 

The correction factors are given in figure 8 for several values of b/c. The dashed 
curves show the theoretical limits expressed in terms of the correction factor F as the 
value of b approaches either zero or infinity. The stress -intensity factor in the 
limiting case (b = °°) was obtained from reference 18 for the case of an edge crack in a 
semi-infinite plate. The other limit is where the elliptical hole reduces to a crack. The 
correction factors approached unity for all finite values of b/c as the crack length 
approached infinity. The edge-crack solution approached 1.12. For b/c ratios greater 
than 2 and large values of a/c the results indicated that the influence of the elliptical , 
hole may not be negligible. However, for small b/c ratios the correction factor 
approached unity very rapidly. 


Crack in a Finite Plate 


In the following section the stress intensity solutions for two cases of a crack in a 
rectangular plate are presented.. Double precision was used in the computer solution of 
the resulting equations primarily because the least-squares procedure generated matrix 
elements whose values ranged over many orders of magnitude and the unknown coeffi- 
cients may have been susceptible to round-off error, hi fact, for some choices of the 
pertinent parameters the single-precision routine did not appear to converge. However, 
for the same parameters the double-precision routine did appear to converge as the 
number of coefficients increased. 


Crack subject ed to wedge-force loading.- For the case of a crack in a rectangular 
plate subjected to wedge-force loading on the crack surfaces (see fig. 9) the stress func- 
tions were taken to be 


Pa 






( 12 ) 


Pa 


fi(z) = Oo(z) + ^ K 

27tzVz^ - 


J 


where 4>q and Oq are given by equations (AlO). The additional terms in equations (12) 
were added in order to account for the concentrated forces on the crack surfaces (ref. 9). 
The unknown coefficients from equations (AlO) in equations (12) were evaluated by 
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satisfying the conditions on the external boundary. To show convergence, an example 

problem with = 0.9 was solved for several values of 2N (see fig. 9). This example 

w 

problem was selected because the close proximity of the crack tip to the boundary was 
expected to pose difficulties in achieving convergence. However, the convergence was 
rapid. 

The correction factors for several values of H/W are shown in figure 10. For a 
given value of H/W the correction factors ranged from unity, at small crack lei^th to 
plate width ratios, to very large values at large ratios. 

In figure 11, a comparison is made between experimental and theoretical stress- 
intensity correction factors for a wedge-force loaded panel. The material was 
7075-T6 aluminum alloy cycled at constant amplitude with a load ratio (minimum load to 
maximum load) of 0.05. The experimental stress-intensity factors were obtained from 
measured crack growth rates used together with a curve of stress intensity as a function 
of crack growth rate (ref. 4). The agreement between the calculated and the experimental 
correction factors is considered good. 

Cracks emanating from a c irctilar hole subjected to uniaxial stress .- For the case 
of cracks growing from a circular hole in a rectangular plate subjected to a uniaxial 
stress (see fig. 12) the stress functions were taken to be the sum of equations (AlO) and 
(All) since both internal and external boundaries must be considered. The unknown coef- 
ficients were evaluated by satisfying the conditions on both boundaries. 

The correction factors for several values of 2R/W and a plate aspect ratio 

(H/W) of 2 are shown in figure 12. As the hole radius approaches zero, the boundary- 

value problem reduces to that of a single crack in a rectangular plate. In this case all 

coefficients in equation (All) are set equal to zero. The correction factors shown 
2R 

for = 0 agree almost identically with Isida's solution (ref. 20) for the case of a 
w 

central crack in a finite width strip. The solid curve was slightly higher than Isida's 

Oo 

solution for values of -^ > 0.7. However, the difference may be due to the influence of 

w 

the finite height. 


CONCLUDING REMARKS 

An improved method of boxmdary collocation was presented and applied to the stress 
analysis of cracked plates having various boundary shapes and subjected to inplane 
loading. The solutions were based on the complex variable method of Muskhelishvili and 
a modification of the numerical method of boundary collocation. The complex series 
stress functions formulated for cracked plates automatically satisfy the boundary condi- 
tions on the crack surfaces. The conditions on the other bovindaries were satisfied 
approximately by the series solution. Stress -intensity correction factors were presented 
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for several configurations involvii^ cracks in the presence of stress concentrations or 
boundaries. The least-squares botindary-collocation method used to minimize the 
resultant -force residuals on the boundary gave better numerical convergence in the 
boundary conditions than the other three collocation methods investigated. The other 
methods were (1) direct specification of boundary forces, (2) direct specification of 
boundary stresses at discrete points, and (3) least-squares specification of boundary 
stresses. The improved method was able to analyze more complex boundary -value 
problems than the other methods because of the more rapid convergence in the series 
solution. 

The stress-intensity factors presented in this paper cover a moderate range of 
configurations and loadings. Most of the botmdary -value problems solved do occur fre- 
quently in aircraft structures such as cracks growing from or near cutouts (windows, 
rivet holes, etc.) or from other structural discontinuities. Thus, the solutions presented 
here may be used to design a structure more efficient against fatigue and fracture or to 
help monitor structural damage due to fatigue loading. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 15, 1971. 


14 



APPENDIX A 


FORMULATION OF THE COMPLEX SERIES STRESS FUNCTIONS 
FOR CRACKED BODIES 


One of the major advances in the field of two-dimensional linear elasticity has been 
the complex-variable approach of Muskhelishvili (ref. 16). The representation of bihar- 
monic functions by analytic functions has led to a general method of solving plane-strain 
and generalized plane-stress problems. 

The following formulation for two-dimensional cracked bodies is based on the work 
of Muskhelishvili for an infinite plane region containing cracks and subjected to inplane 
loading (see fig. 13). The known surface tractions are applied to the boundaries of this 
region. The body forces are assumed to be zero and the material is assumed to be linear 
elastic, isotropic, and homogeneous. The equilibrium and compatibility equations are 
combined to form the biharmonic equation 


3^U , o a^u ^ e^u _ Q 

9x2 ay2 8y4 


(Al) 


where U(x,y) is the Airy stress function. 

The biharmonic function U(x,y) can be ejq>ressed as 


U(x,y) = ^ Re y 'F(z)dz + y $(z)dz + (z - z)^(z) 


(A2) 


where #’(z) = <p{z) and ^’(z) = S2(z). The primes and bars denote differentiation and 
complex conjugates, respectively. Therefore, the generalized plane-stress and plane- 
strain problems are reduced to the determination of these functions subject to specified 
boundary conditions. 

From reference 16, the stress functions in the neighborhood of a crack (z = ±a) 
can be expressed as 


^(z) 

fi(z) 


Ql(z) 

^z2 - a2 


+ Q2(z) 


Q0(z) 

\/z2 - a2 


+ Q4(z) 




(A3) 
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APPENDIX A - Continued 


where 

oo 

Qs(^) “ ^ ■^sn^^ (■A.4) 

n=-co 

In the present investigation, the formulation of the complex series stress functions 
for simply and multiply connected regions containing cracks was restricted to the situa- 
tion where the configuration and loading are symmetric about the x- and y-axis. The 
boundary conditions satisfied by (p S2 for a stress-free crack surface are 

(I) ^y = '^xy = 0 (|x|< 

(II) 'Txy = v = 0 (|x|g 

(III) Txy=u = 0 (|y|^ 

where the coordinate system used is shown in figure 13. 
of the stress functions are 

+ ay = 2 [^(z) + 0(z)J 

ay - o-x + 2iTxy = 2 [(z - z)0'(z) - 4>{z) + Jl(z^ 

The displacements are given by equation (1). 

The boundary conditions in equations (A5) define a unique relationship between the 
two analytic functions (p and The relationship for the series which contains the 
square-root term is QjCz) = Q 3 (z) and for the series which contains no square-root 
term, Q 2 (z) = -Q 4 (z)- These conditions insure the satisfaction of the boundary conditions 
(stress free) on the crack surfaces. 

For a multiply connected region, in addition to the boundary conditions stated in 
equations (A5), the single-valuedness of displacements must be insured. This condition 
is stated as 


a; 

II 

a; 

II 

0; 

"o" 

11 


J 


(A5) 


The stresses expressed in terms 


(A6) 


K ® (i(z)dz - ® Sl(z)dz = 0 

Ja. 


(A7) 


where Aj is the contour around each separate hole boundary Lj (see fig. 13). 
The stress intensity factor K is obtained from the stress functions as 


16 



APPENDIX A - Continued 


K = 2 lim \jz - a cp{z) 
z— a 


(A8) 


As an example, consider a crack located along the x-axis in an infinite plate as 
shown in figure 14. The dashed lines Lq and Lj define the boundaries of the region 
(shaded area). The internal boundary has cracks growing from the edge of the hole 
to X = ±a. The boundaries Lq and Lj may have any simple shape which is ssmnmetric 
about the x- and y-axis and be subjected to any symmetric boundary conditions. The 
complex-series stress functions must be single-valued and analytic in the region between 
Lq and Lj. This region does not include the crack which is represented by a branch 
cut. The stress functions are 


(p{z) = (^q(z) + (piiz) 
S2(z) = ^2 q(z) + Oi(z) 


(A9) 


The subscripts denote the functions which are used to satisfy the boundary conditions on 
boundaries Lq and Lj, respectively. The stress functions used to satisfy the condi- 
tions on the external boundary Lq are 




N N 

n=0 n=0 


(AlO) 


where the 
is located 
are zero. 


coefficients 
at infinity 


Bji are real. 


Ajj and 

Ao=f, Bo='f(X-l) 


In the situation where the boundary Lq 


2’ 4 

For the internal boundary Lj 


and the remaining coefficients for n = 1 
the stress functions are 




N 




n=l 



1=1 



(All) 


where the coefficients Cj^ and are real. The stress functions, equations (AlO) and 
(All), automatically satisfy the boundary conditions on the crack surfaces. The conditions 
on the other boundaries (Lq and Lj) were approximated by the series solution. The 
stress-intensity factor K calculated from equations (A8), (AlO), and (All) is 


K = /?a 



2Aj^a2n + 


N 


Y 

Z/ n2n 


n=l 



(A12) 
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APPENDIX A - Concluded 


In the case of multiple circular holes or elliptical holes, the complex- series stress 
functions require the use of poles at various stations along the x- or y-axis. The stress 
functions are given by 



(A13) 


where the coefficients and Djn are real. In these stress functions the poles ±Zj 

must lie on either the x- or y-axis and be symmetrically placed about the other axis. 
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APPENDIX B 


OTHER COLLOCATION METHODS INVESTIGATED 


Direct Specification of Boundary Stresses 


The boundary- collocation method treated in this section involves the specification of 
the normal- and shear-stress components at discrete points on the boundary. The 
boundary- value problem considered had cracks emanating from the edges of a circular 
hole in an infinite plate subjected to a biaxial stress (see fig. 15). The complex equation 
for the two stress components on the boundary is 


% “ i'^nt = + (p(z) - [(z - z)0’(z) - 4>{z) + 

where the stress functions are 




I 


n=0 


SiL 

z2n 


± 


N 


I 

n=0 



(Bl) 


(B2) 


The coefficients Cq and Dq were determined from the stress conditions at infinity. 
The remaining coefficients were determined from the conditions that and = 0 

at N equally spaced points on the boundary. The resulting equations were solved on a 
computer using single precision. The stress-intensity factor was calculated from 
equation (A8) as 


K 



n=l 




(B3) 


where the term in the brackets is the correction factor F, 


Least-Squares Specification of Boundary Stresses 

The boundary-collocation method described in this section is similar to the one in 
the preceding section, except that here more equations are written in terms of the stress 
boundary conditions than there are unknown coefficients. The coefficients are then chosen 
so as to minimize the sum of the stress residuals on the boundary at a specified number 
of points M. The normal stress and shear stress components expressed in terms of the 
stress functions (eqs. (B2)) are 
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APPENDIX B - Concluded 


N N 

“ X ^ ^n^2n 

n=0 n=0 


'nt 


N N 

“ ^ ^n^ln + ^ ^nG2n 


n=0 


n=0 


(B4) 


where Fg^j and Gs^ (s = 1,2) are the influence functions derived from the stress func- 
tions for the respective coefficient. Again, the coefficients Cq and Dq are deter- 
mined from the stress conditions at infinity. Since the circular boundary is stress free, 
the square of the error in the boundary condition at point is 


'm 


fv 

N 


"r 

N 1 

A C„Fl„ + 

2y ^n^2n 


^ CnGiji + 

) DnG2n> 

n=0 

n=0 

J 

1 

n=0 

n=0 J 


m 


m 


(B5) 


The coefficients were then selected so that the sum of the squares of the errors at a num 
ber of points M on the boundary was a minimum 


M 

= 0 

9Cq 

M 

=0 

J 


(B6) 


where q = 1,2,. . . ,N. Equation (B6) results in a set of 2N linear algebraic equations 
for Cj^ and Dj^ 
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(B7) 


N 

N 


2 j “nq^n 

1 

^nq^n ~ ® 

n=0 

n=0 

> 

N 

N 


X ^qnCn 

t 

>o 

11 

O 

n=0 

n=0 

J 


where 


O' 


nq - 


M 

^ ^In^lq + Gln^lq) 




nq 


M 


^2n^lq G2nGiq^^ 


’'nq Z ^2n^2q ^2n‘^2q) 

m=l ^ 

These equations were solved on a computer using single precision. The resulting coeffi- 
cients were used in equation (B3) to calculate the stress-intensity factor. 


21 



REFERENCES 


1. Paris, Paul C. ,* Gomez, Mario D.; and Anderson, William E.: A Rational Analytic 

Theory of Fatigue. Trend Eng. (Univ. of Washington), vol. 13, no. 1, Jan. 1961, 
pp. 9-14. 

2. Donaldson, D. R. ; and Anderson, W. E.: Crack Propagation Behavior of Some Airframe 

Materials. D6-7888, Paper presented at Crack Propagation Symposium 
(Cranfield, England), Sept. 1961. 

3. Anon.: Fracture Testing of High-Strength Sheet Materials: A Report of a Special 

ASTM Committee. ASTM Bull., no. 243, Jan 1960, pp. 29-40. 

4. Figge, I. E.; and Newman, J. C. , Jr.: Prediction of Fatigue-Crack-Propagation 

Behavior in Panels With Simulated Rivet Forces. NASA TN D-4702, 1968. 

5. Bowie, O. L. : Analysis of an Infinite Plate Containing Radial Cracks Originating at the 

Boundary of an Internal Circular Hole. J. Math. Phys., vol. XXXV, no. 1, Apr. 1956, 
pp. 60-71. 

6. Grebenkin, G. G.; and Kaminskii, A. A.: Propagation of Cracks Near an Arbitrary 

Curvilinear Hole. Sov. Mater. Sci., vol. 3, no. 4, July-Aug. 1967, pp. 340-344. 

7. Fichter, W. B. : Stresses at the Tip of a Longitudinal Crack in a Plate Strip. NASA 

TR R-265, 1967. 

8. Isida, M.: On the Determination of Stress Intensity Factors for Some Common 

Structural Problems. Eng. Fracture Mech. , vol. 2, no. 1, July 1970, pp. 61-79. 

9. Erdogan, Fazil: On the Stress Distribution in Plates With Collinear Cuts Under Arbi- 

trary Loads. Proceedings of the Fourth U.S. National Congress of Applied 
Mechanics, Vol. One, Amer. Soc. Mech. Eng., 1962, pp. 547-553. 

10. Hulbert, L. E.; Hahn, G. T.; Rosenfield, A. R.; and Kanninen, M. F.: An Elastic- 

Plastic Analysis of a Crack in a Plate of Finite Size. Applied Mechanics, 

M. Het^nyi and W. G. Vincenti, eds.. Springer- Verlag, 1969, pp. 221-235. 

11. Kobayashi, A. S.; Cherepy, R. D. ; and Kinsel, W. C.: A Numerical Procedure for 

Estimating the Stress Intensity Factor of a Crack in a Finite Plate. Trans. ASME, 
Ser. D.: J. Basic Eng., vol. 86, no. 4, Dec. 1964, pp. 681-684. 

12. Bowie, O. L. ; and Neal, D. M. : A Modified Mapping- Collocation Technique for Accurate 

Calculation of Stress Intensity Factors. Int. J. Fracture Mech., vol. 6, no. 2, 

June 1970, pp. 199-206. 

13. Conway, H. D. : The Approximate Analysis of Certain Boundary- Value Problems. 

Trans. ASME, Ser. E: J. Appl. Mech., vol. 27, no. 2, June 1960, pp. 275-277. 


22 



14. Hulbert, Lewis Eugene: The Numerical Solution of Two-Dimensional Problems of the 

Theory of Elasticity. Bull. 198, Eng. Exp. Sta., Ohio State Univ. 

15. Hooke, C. J.: Numerical Solution of Plane Elastostatic Problems By Point Matching. 

J. Strain Anal., vol. 3, no. 2, Apr. 1968, pp. 109-114. 

16. Muskhelishvili, N. I. (J. R. M, Radok, transl.): Some Basic Problems of the Mathemat- 

ical Theory of Elasticity. Third ed., P. Noordhoff, Ltd. (Groningen), 1953. 

17. Fox, L.: An Introduction to Numerical Linear Algebra. Oxford Univ. Press, 1965. 

18. Paris, Paul C.; and Sih, George C.: Stress Analysis of Cracks. Fracture Toughness 

Testing and Its Applications, Spec. Tech. Publ. No. 381, Amer. Soc. Testing 
Mater., c.1965, pp. 30-83. 

19. Savin, G. N. (Eugene Gros, transl.): Stress Concentration Around Holes. Pergamon 

Press, Inc., 1961. 

20. Isida, M.: On the Tension of a Strip With a Central Elliptical Hole. Trans. Japan Soc. 

Mech. Engrs., vol. 21, no. 107, 1955, p. 511. 


23 


TABLE I.- CRACKS EMANATING FROM A CIRCULAR HOLE IN AN 
INFINITE PLATE SUBJECTED TO BIAXIAL STRESS 



TABLE n.- CRACKS EMANATING FROM AN INTERNALLY 
PRESSURIZED HOLE IN AN INFINITE PLATE 



a/R 

1.01 

1.02 

1.04 

1.06 

1.08 

1.10 

1.15 

1.20 

1.25 

1.30 

1.40 

1.50 
1.60 
1.80 
2.00 
2.20 

2.50 

3.00 

4.00 


Fl (A=l) 

Fq (A=0) 

0.2188 

0.1725 

.3058 

.2319 

.4183 

.3334 

.4958 

.3979 

.5551 

.4485 

.6025 

.4897 

.6898 

.5688 

.7494 

.6262 

.7929 

.6701 

.8259 

.7053 

.8723 

.7585 

.9029 

.7971 

.9242 

.8264 

.9513 

i .8677 

.9670 

.8957 

.9768 

.9154 

.9855 

.9358 

.9927 

.9566 

.9976 

.9764 



TABLE ni.- CRACK LOCATED BETWEEN TWO CIRCULAR HOLES 
m AN INFINITE PLATE SUBJECTED TO UNIAXIAL STRESS 


s 



K = S F 



^ (1=6) 
0,8661 
.8671 
.8700 
.8748 
.8810 
.8969 
.9150 
.9328 
.9486 


F (|=10) 
0.9506 


.9512 


.9528 

.9553 

.9585 

.9622 

.9663 


26 



TABLE IV.- CRACK APPROACHING TWO CIRCULAR HOLES IN AN 
INFINITE PLATE SUBJECTED TO UNIAXIAL STRESS 
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Figure 8.- Stress -intensity correction factors for cracks emanating from 
an elliptical hole in an infinite plate subjected to uniaxial stress. 



















